Typeset with jpsj2.cls <ver.l.2> 



Full Paper 



Temperature dependence of ESR intensity for the nanoscale molecular 

magnet V15 

Manabu Machida^ *, Toshiaki Iitaka^ \ and Seiji MlYASHlTA^ * 

^Institute of Industrial Science, The University of Tokyo, 4-6-1 Komaba, Meguro-ku, Tokyo 

153-8505, Japan 

'^Computational Astrophysics Laboratory, RIKEN (The Institute of Physical and Chemical Research), 

2-1 Hirosawa, Wako, Saitama 351-0198, Japan 
^Department of Physics, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, 

Bunkyo-ku, Tokyo 113-0033, Japan 
(Received February 2, 2008) 

The electron spin resonance (ESR) of nanoscale molecular magnet V15 is studied. Since the 
Hamiltonian of V15 has a large Hilbert space and numerical calculations of the ESR signal 
evaluating the Kubo formula with exact diagonalization method is difficult, we implement 
the formula with the help of the random vector technique and the Chebyshev polynominal 
expansion, which we name the double Chebyshev expansion method. We calculate the tem- 
perature dependence of the ESR intensity of V15 and compare it with the data obtained in 
experiment. As another complementary approach, we also implement the Kubo formula with 
the subspace iteration method taking only important low- lying states into account. We study 
the ESR absorption curve below lOOK by means of both methods. We find that side peaks 
appear due to the Dzyaloshinsky-Moriya interaction and these peaks grows as temperature 
decreases. 
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The V15 molecule is the complex of formula Kg [V^yAs6042 (H2O)] • 8H2O. Since Miiller 
and Doring synthesized this molecule for the first time,^"^ V15 has been studied intensively as 
one of promising nanometer-scale molecular magnets. In V15, fifteen 1/2 spins of vanadium 
ions are arranged almost on a sphere. The triangle cluster with three spins is sandwiched by 
the upper and lower hexagons. In experiment, an adiabatic change of the magnetization has 
been observed in a fast sweeping field. In a slow sweeping field, an interesting magnetic 
plateau appears. This phenomenon, called the phonon bottleneck effect, is due to an effect of 
contact with thermal bath,^' and also theoretically analyzed from a general viewpoint of the 
magnetic Foehn effect.^ Since the magnetization changes smoothly at zero field from —1/2 to 
1/2 when the field is swept adiabatically at low temperatures, the gap between the ground 
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state and the lowest excited state at zero field is believed to open. The most plausible origin of 
the gap is the Dzyaloshinsky-Moriya (DM) interaction.^' ^'^"^^ Around 2.8T, the ground state 
magnetization changes from 1/2 to 3/2. This change also occurs smoothly, and this broadness 
of the change is also considered to be caused by the DM interaction.'' However, the detail of 
the mechanism of this broad change is not fully understood yet.^^'-^^ 

One possible way to determine the mechanism of adiabatic change in V15 is that we com- 
pare the numerically obtained ESR absorption curve of V15 with the experimentally obtained 
one. With this motivation, in this paper, we establish numerical methods to study the ESR 
absorption curve obtained from the model Hamiltonian of V15. The ESR absorption curve is 
calculated by the Kubo formula in theory.^ The direct implementation of the Kubo formula 
is to diagonalize the Hamiltonian matrix. However, this direct method is impossible in most 
cases because the dimension of the Hilbert space is quite large. For example, it is 2^'^ = 32768 
for V15. 

In this paper, we study the ESR absorption curve of V15 by proposing two methods: one is 
the double Chebyshev expansion method (DCEM) and the other is the subspace iteration method 
(SIM). Relying on these methods, we reveal the temperature dependence of the intensity of 
V15, and find side peaks due to the DM interaction which allows transitions between excited 
states otherwise forbidden. 

The magnetic interactions in V15 are described by the following Hamiltonian^^' 

(id) * 

We show the interactions between spins in Fig. 1. For Jij, we have three different values J, 
Ji, and J2 (|J| > IJ2I > 1^1 1) with respect to the bonds on the upper and lower hexagons. 
Three spins on the triangle between two hexagons interact with the hexagons by Ji and J2. 
The interactions between the three spins are negligibly small. Here we take J = —BOOK, 
J2 = — 350K, and Ji = — 225K.^^ The second term on the right-hand side in Eq. (1) describes 
the DM interaction. DM vectors are considered to exist on the two hexagons in the bonds 
with J. We take the reference DM vector Di^2 to be Df ^ = Df ^ = ^!,2 = The 
other DM vectors on the upper hexagon are obtained by rotating -Di,2 by 27r/3 and 47r/3, 
i.e., 4 = 14.641K, = -54.641K, = 40K, Df g = -54.641K, L>|g = 14.641K, and 
D^Q = 40K. If we assume the symmetry of Vis,^ the lower hexagon differs from the upper 
hexagon by rotation 7r/6, and the z components of the DM vectors on the lower hexagon have 
opposite sign from those of the DM vectors on the upper hexagon. Thus, the DM vectors on 
the lower hexagon arc obtained as Df^ = — 54.641K, D^q = — 14.641K, D^q ^-^ = — 40K, 
^12,13 = 40K, 1)^2,13 = -40K, Z)f2,i3 = -40K, L>f4 ^5 = 14.641K, Df^^^r^ = 54.641K, and 
Df^ 15 = — 40K. We assume the magnetic field is applied parallel to the c-axis of the molecule 
(z-axis) and set Hs = (0, 0, Hs). We take Hs = 4T throughout the paper (see also Fig. 4(b)). 
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Fig. 1. The figures on the left show interactions within each hexagon. The interaction within the 
triangle is negligibly small. The figure on the top right shows interactions between the triangle 
and the upper and lower hexagons. 

Double Chebyshev expansion method (DCEM). — This method realizes 0{N) calculation 
with the help of the random vector technique^'^''^^ and the Chebyshev polynomial expansion of 
exponential operators. The DCEM is an extension of the Boltzmann-weighted time-dependent 
method (BWTDM) developed by litaka and Ebisuzaki.^^ We briefly review the BWTDM. 
The procedure of the BWTDM is divided by the following five steps. By adopting the unit in 
Ref. 22, we have dimensionless energy, time, and temperature. At the first step, we prepare a 
random vector |<^). For a given basis \n) of the Hilbert space, this random vector is given by 
|<I>) = Yln=i 1^) ^n- Here, the dimension of the Hilbert space is N and the statistical average of 
the complex random coefficients {^„} satisfies = (5n'n- At the second step, we obtain 

the Boltzmann-weighted random vector |$Boitz) = e~^^/^ |<I>) making use of the Chebyshev 
polynomial expansion. Here, /?(= 1/T) is the inverse temperature. To this end, we divide the 
Hamiltonian by constant A A as TYsc = Ti/AX so that the largest eigenvalue of Tisc does not 
exceed unity. Then, we expand e~^^/^ as follows. 

e-'^w/^ = /o(-/3AA/2)ro(Hsc) 

+ 2^4(-/3AA/2)rfe(?isc), (2) 

k=l 

where Ikix) is the modified Bessel function and Ti^iTisc) is the Chebyshev polynomial, which 
satisfies Tk(Hsc) = 2HscTk~i{'Hsc) - Tk^2{'Hsc), To(Hsc) = 1, and Ti{7isc) = Hsc- By this 
procedure, we obtained e~^^/^ without diagonalization. At the third step, we obtain time 
evolutions of vectors |$Boitz;i) = e~'^* |$Boitz) and \^M^',t) = e~''^* |$Afa:), where \^m^) = 
I^Boitz) and = Ylj'=T ^j- ^'^ BWTDM, the time evolution is performed by the 
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leap frog method, which evolves state \(t>\t) as 

t + At) = -2inAt 10; t) + !</.; t - At) . (3) 

Note that the condition E^maxAt ^ 1 should be satisfied, where E^ax is the largest eigenvalue 
of the Hamiltonian. At the fourth step, we calculate the correlation function 



git;T) = {M^M^t)) 



= TV 



/TV 



(( ($M^;t|M^|$Boitz;t) )) . . 

(( ($Boltz|$Boltz) )) ' ^ ' 

where the trace was replaced by the inner product of random vectors. Finally, the imaginary 
part of the dynamical susceptibility x"i^jT) is obtained by the Fourier transform of g{t;T). 

X"(a;; T) = (l - e"^-) Re g{t; T)e~'^'dt 

l-e-^'^jRey 5(i;T)e-''"*e-'' * /^di. (5) 

Here, we introduced the Gaussian filter with variance 1/??^. This r] determines the frequency 
resolution. The upper limit of the integral Tmax satisfies Tmax ~ 1/?? in order to avoid the 
Gibbs oscillation. Also r] should satisfy < ^ 1, <C Ha,, and firf -C Ha,. The average 
energy absorption per unit time I{u);T) is given by 

7(a;;r) = ^^X>;r). (6) 

The DCEM is almost the same as the BWTDM. Only the third step is different. In the 
DCEM, we make use of the Chebyshev polynomial expansion not only in the second step 
obtaining o'^'^l'^ , but also in the third step of time evolution. State vector \4>\t) is evolved 
with the Chebyshev polynomials as 

|<^;i + r) = e--^^«-|<^;i) 

= Jo(rAA)ro(Hsc)|<^;t) 

+ 2 ^ (-if Jfc(TAA)Tfc(?^sc) (7) 
fe=i 

where {x) is the Bessel function. Note that time step r is not necessarily small. In the ESR 
experiment for V15, magnetic field Ha, (~ IK) is usually much smaller than the strongest 
coupling Jmax(= \ J\) between spins (~ lO^K). That is why the frequency of precession of the 
spins is rather small. This means that we need to evolve state vectors long time but do not 
need fine resolution of time step in order to detect small frequencies of spin precession. Since it 
is possible for time step r to have larger value than time step At, the DCEM is more efficient 
than the BWTDM for low magnetic fields. Table I shows typical computation times with the 
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Table I. Typical computation times of two methods at various magnetic fields. 
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Fig. 2. The temperature dependence of the intensity of V15 is shown. Solid circles denote data by 
the DCEM and open squares denote data by the SIM with the lowest eight levels. Dashed lines 
denote, from the bottom, intensities of independent 1, 2, 3, and 15 spins, respectively. 



DCEM and the BWTDM. Using asymptotic behavior of the Bessel function, the ratio of two 
computation times is evaluated as 

[Chebyshev] oi/s 
[Leap-frog] Jmax 

where a and h are constants. 

We first study the temperature dependence of the intensity I{T) (= I{oj;T)<1oj) us- 
ing the DCEM. In Fig. 2, we show the temperature dependence of I{T). The intensity I{T) 
obtained by the DCEM is denoted by solid circles. The open squares in Fig. 2 denote data 
obtained by the subspace iteration method, which we will explain below. At very high tem- 
peratures, each spin in V15 acts as an isolated spin. Therefore, the intensity should be given 
by the multiplication of Ii{T) by 15, where Ii{T) = ^H'^Hq, ianh. {(3 H'^/ 2) is the intensity of 
an isolated spin. In Fig. 2, 15 x Ii{T), 3 x /i(T), 2 x Ii{T), and 1 x Ii{T) are also shown 
by dash-dotted line, dotted line, short-dashed line, and dashed line, respectively. As temper- 
ature decreases, the effective number of spins changes from 15 to 3. Note that the ground 
state magnetization is 3/2 in the present case of = 4T (see also the energy diagram in 
Fig. 4(b)). In experiment, the intensity starts deviating the curve for three spins at around 
200K.24 Here we confirm that the intensity actually follows the curve for independent 15 spins 
at high temperatures. 

Subspace iteration method (SIM). — We explain another method obtaining the ESR ab- 



Hs/Jn 



(8) 
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sorption curve of large molecules. At low temperatures, the ESR absorption occurs only from 
transitions among low-lying energy levels because states of the system are confined near the 
ground state. The SIM is a numerical method that takes a small subspace of the total Hilbert 
space which concerns only low-lying states. ^^'^^ Within this subspace, the ESR absorption 
curve is obtained by the explicit formulation of the Kubo formula. That is, by direct diagonal- 
ization of the small Hamiltonian. We obtain the eigenvalues and eigenvectors of the subspace 
and obtain the imaginary part of the susceptibility x"i^)i^^ 

m,n 

xS{u;-{En-Em)), (9) 

where Z is the partition function. 

The calculation in the DCEM is numerically exact, but it generally requires longer com- 
putation time than that of the SIM. On the other hand, the calculation in the SIM is carried 
out shorter time for a moderate value of N, which is enough to study the property at low tem- 
peratures. However, it ignores higher states completely. The two methods are complementary 
to each other. 

In Fig. 2, the intensities obtained by the SIM with the lowest eight levels (N = 8) are 
plotted by open squares. Those are consistent with the solid circles obtained by the DCEM 
at temperatures lower than 200K. Although the data of the DCEM start deviating from the 
line of three spins at around 200K, the data of the SIM stay at the line, which reflects the 
fact that the SIM considers only the lowest eight levels. Here we confirm that the SIM can 
reproduce the data of the DCEM at low temperatures with a shorter computation time. Thus, 
in Fig. 2, we can have more points at low temperatures. 

We point out the following two points, which will be discussed in more detail elsewhere. 
First, the data slightly deviate from the dotted line at low temperatures as seen in the exper- 
iment by Sakon et al.^^ For the present field {Hs = 4T), i.e., the ground state magnetization 
is 3/2, the ratio /(T)//i(T) takes values smaller than 3 at low temperatures due to the DM 
interaction. Second, in the experiment by Ajiro, et al.,^^ the intensity follows the line of an 
isolated 1/2 spin {I{T)/Ii{T) = 1). This is because the ground state magnetization is 1/2 in 
the field weaker than 2.8T. 

Finally, let us study the absorption curve I{u);T). Figure 3 compares I{uj;T = 32K) and 
I{io; T = 64K) obtained by the DCEM. Note that the widths of the peaks are solely due to the 
finite frequency resolution of the DCEM given by r]. We see that the peaks grow as temperature 
decreases. In the upper panel in Fig. 3, intensities with and without the DM interaction are 
compared. We conclude that the side peaks near the main peak around 110GHz(~ 4T) appear 
due to the DM interaction. 

Figures 4(a) and 4(b) show results obtained by the SIM. In Fig. 4(a), the ESR absorptions 
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Fig. 3. The upper and lower panels show the ESR absorption curves of V15 at 4T and at 32K and 
64K, respectively. In the upper panel, the dashed curve shows the absorption curve obtained from 
the Hamiltonian without the DM interaction. 




Fig. 4. (a) The ESR absorptions calculated by the SIM at 4T and 32K are shown. The height of 
the pulse at each resonant frequency expresses the value of the coefficient of the delta function of 
the absorption. The inset shows the magnified figure around llOGHz. (b) The lowest eight energy 
levels of V15. The vertical dashed line at Hs — 4T shows the magnetic field applied to V15. The 
arrows denote transitions corresponding peaks labeled by a through e in (a). 



by the lowest eight levels (A^ = 8) are calculated at 32K. The height of pulses shows the value 
of coefficients of delta functions of the absorptions. The peak around llOGHz consists of three 
resonant peaks, whereas each side peak consists of single resonant peak. Figure 4(b) shows 
the lowest eight energy levels of V15 as a function of field Hs- Each labeled peak in Fig. 4(a) 
comes from the corresponding labeled transition between energy levels shown in Fig. 4(b). 

Thus, the combination of the DCEM and the SIM provides a powerful method to elucidate 
the ESR of large molecular magnets such as V15. 

This work is supported by the Grant-in-Aid from the Ministry of Education, Culture, 
Sports, Science and Technology, and also by NAREGI Nanoscience Project, Ministry of Edu- 
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